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It is shown that thermal fluctuations of the normal component induce dephasing - reversible damp- 
ing of the low energy collective modes of a confined Bose-Einstein condensate. The dephasing rate 
is calculated for the isotropic oscillator trap, where Landau damping is expected to be suppressed. 
This rate is characterized by a steep temperature dependence, and it is weakly amplitude dependent. 
In the limit of large numbers of bosons forming the condensate, the rate approaches zero. However, 
for the numbers employed by the JILA group, the calculated value of the rate is close to the exper- 
imental one. We suggest that a reversible nature of the damping caused by the thermal dephasing 
in the isotropic trap can be tested by the echo effect. A reversible nature of Landau damping is also 
discussed, and a possibility of observing the echo effect in an anisotropic trap is considered as well. 
The parameters of the echo are calculated in the weak echo limit for the isotropic trap. Results of 
the numerical simulations of the echo are also presented. 



o 



o 
o 



> 

00 



I 

O 

o 



X 



PACS numbers; 03.75.Fi, 05.30.Jp, 32.80.Pj, 67.90.-(-z 



I. INTRODUCTION 

Properties of the trapped ultracold atomic gases 
demonstrating the phenomenon of Bose-Einstein conden- 
sation ||l|-|[ attract much of attention. A confined geom- 
etry of the cloud leads to circumstances under which new 
phenomena such as, e.g^ a specific scaling in the thermo- 
dynamical properties [^ and the two step condensation 
g can be anticipated. It was also suggested that the 
phase of the confined condensate should demonstrate the 
phenomenon of collapses and revivals pi and the phase 
diffusion effect |j]. 

Recently it has been discussed in Refs. MM, that the 
dynamical quantum behavior of the atomic cloud con- 
sisting of a finite number of bosons is very different from 
the properties of infinite systems. Specifically, an ampli- 
tude of a normal mode of the confined condensate should 
demonstrate quantum dephasing which results in an ap- 
parent damping of the mode even at zero temperature T. 
The rate of such damping is determined by the interpar- 
ticle interaction and turns out to be linearly dependent 
on the mode amplitude, if it is small MM- A possibility 
of the mode revival has also been predicted ^^ . A cur- 
rent experimental situation does not allow to resolve this 
effect because of the presence of the normal component 
which introduces a substantial damping on its own. 

Damping of the normal modes at finite T of the con- 
fined condensate has been experimentally studied by the 
JILA and the MIT [|ll| groups. Recently it was 
stressed by Pitaevskii ||l5] (see also in Refs. O-Ill]) that 
the damping in the trap containing condensate is essen- 
tially Landau damping (LD) which is coUisionless in na- 
ture. In other words, it is rather the dephasing of the 
collective modes than their irreversible dissipation [[l7| . 
It is worth noting that the reversible nature of the LD 
in the classical uniform plasma can be revealed in the 
plasma echo effect 11^ . 

Theoretical approaches |l^-|lq| are based on applying 



the standard results of the perturbation theory developed 
for an infinite and uniform medium. Accordingly, a dis- 
crete structure of the quasiparticle spectrum is tacitly 
replaced by an effective continuum |[l^-[lq| . In the recent 
work |l^ it has been emphasized that the existence of the 
LD is directly related to a presence of the randomness in 
the spectrum of the anisotropic traps. Conversely, in the 
isotropic trap the LD should be suppressed because the 
quasiparticle spectrum is regular [ [l8[ . 

In this paper, we suggest that in a confined condensate 
where the LD is suppressed the damping of the collective 
modes can still be observed. Specifically, thermal fiuc- 
tuations of the population factors of the normal compo- 
nent induce a reversible dephasing of the collective os- 
cillations. The nature of such a damping is similar to 
that of the dephasing observed in quantum dots (see in 
Ref. |jl9]). We show that, while being essentially zero 
for a traditional uniform condensate, the rate of such a 
dephasing I/t^ in the trap containing 10^ — 10^ atoms 
can be comparable to the rate of the damping observed 
experimentally in Ref. pO| . 

We also suggest that reversibility of the damping in a 
confined condensate can be tested in a kind of the phonon 
echo experiment (see in Ref. pO[), when two consecutive 
external pulses imposed on the trap induce a third pulse 
- the echo - at the time approximately equal to twice 
the time interval between the first two pulses. In this 
paper we analyze the case of the isotropic trap, where 
the LD is not expected to be an important mechanism of 
dissipation. In a future work we will investigate the echo 
in strongly anisotropic traps, where the main mechanism 
of damping is the LD. 



II. THERMAL FLUCTUATIONS AND 
DEPHASING IN THE ATOMIC TRAP 



In the Refs. ISlBl] it has been shown that a collective 
mode of the confined condensate should exhibit a de- 
phasing caused by the inter-particle interaction. In the 
case T = 0, the dephasing is produced by the particles 
forming the mode. Below we will consider a similar effect 
caused by the interaction between a low energy collective 
mode and the quasiparticles forming a normal cloud. 

Let us discuss general reasons for the thermally in- 
duced dephasing in the presence of the inter-particle in- 
teraction. Especially we will clarify why this effect is 
not significant for an infinite condensate, where the only 
cause of the dephasing should be the LD. The following 
analysis is based on the approach suggested by Pitaevskii 
in Ref. Q for the case T = 0. We extend this analysis 
for T^O. 

If an external modulation at some frequency ujq has 
excited a system, the many body time dependent wave 
function \t) constructed in terms of the quasiparticle 
states \No,Ni...) acquires the form 



\t) 



E 

No.Ni. 



CjVo,JVi...e-*^"o. 



'\No,Ni,. 



(1) 



where iV„ is a population number for the n-th {n = 
0, 1, 2, ...) quasiparticle state; Enq^Ni... stands for the en- 
ergy of the state \No, Ni...) and Cno,Ni... denote the nor- 
malized coefficients (J^n n \Cno,Ni...\'^ = !)■ For the 
case of weak interaction between quasiparticles, the ex- 
ternal drive, which is in resonance with the energy of, 
e.g., the 0-th quasiparticle state, creates a coherent mix- 
ture of the quasiparticles in the 0-th state and does not 
affect significantly the other states. Accordingly, one can 
factorize Cno.Ni... as 



C 



No.Ni. 



'No' 



.(1), 



.(2) 

'Nn- 



and assume that 



\C 



(0)|2 



N, 



No 



No\ 



exp(-7Vo), 



(2) 



(3) 



with the rest of the C- coefficients corresponding to the 
number states {C}^ — 1 for some Nn — N^ and Cj^ ~ 
otherwise). In Eq.(3) Nq is given by the value of the 
classical amplitude of the resonant mode SbI . 

The expectation value of the single particle operator 
A{t), which changes the number of quasiparticles in the 
given state by 1, acquires a resonant contribution from 
the 0-th state. This is 



(0)* ^(0) 



{t\A\t) = Ea.„ CZ+iC';;!AeM^{ENo 



+ l,Ni,.. 



Eno.Ni,---)^] + c-c. 



(4) 



where ^ is a corresponding single-particle matrix ele- 
ment. In the case of a weak interaction between quasi- 
particles, one can expand the energy in accordance with 
the Landau theory of quantum liquids as 



Eno,Ni, 



y^ (^nN„ 



^E 



(5) 



where the coefficients gmn arise due to interaction be- 
tween quasiparticles. 

At finite T, the solution (4) should be averaged over the 
initial values of Nn with n > 0, distributed in accordance 
with a thermal ensemble. Substituting (5) into (4) and 
neglecting the term ^ Nq significant for very small T 
only (see [pP]), one obtains 



(<|^|t)^ = Ae-»*^<);i<) 



No 



(exp(it ^ gonN„))T + c.c, 



(6) 



n>0 



with (...)t denoting the thermal average with respect to 
the population factors iV„. 

In what follows we will assume that Nn are distributed 
in accordance with the grand canonical distribution. As 
will be discussed below, this assumption is valid as long as 
the mean number of atoms in the condensate Nc is suf- 
ficiently large. Thus, after straightforward calculations 
one finds 



n>0 



{exp{it'^ gonNn))T = 

1 ~ exp(z.go„t - Un/T) 



(7) 



exp - ^ In -^ 



n>Q 



exp(-w„/T) 



We note that in a large system the matrix elements 
gon are scaled as gon ~ 1/V by the effective volume V of 
the system. Accordingly, in the formal limit 1^ ^ oo, one 
should expand exp(i(7oni) in Eq.(7) up to the first term ~ 
1/V only. Then a smallness of 1/V will be compensated 
by the summation '^^ ~ V. This results in 



(exp(ji ^ gonNn))T = exp(it ^ gonNn), 



(8) 



n>0 



n>0 



where iV„ denotes a Bose distribution of the non- 
interacting quasiparticles. As one can see, in the infinite 
system fluctuations of the numbers of the quasiparticles 
do not cause any dephasing. Instead, the frequency loq 
experiences a thermal shift luq -^ ujq + J2n9onNm and 
the only cause of the dephasing is the LD. 

In the case of a finite system one should keep higher 
order terms in Eq.(7). As will be shown below, the 
dephasing rate in the oscillator trap contains a small- 
ness a/rtrap << 1, where a and rtrap stand for the 



scattering length ^ 10"'^ — 10~^cm and some typical 
scale ^ 10~^cm associated with the trapping poten- 
tial, respectively. Therefore, as long as a/rtrap << 1, 
one can neglect the higher terms o{a'^ /r^^ ). Such an 
approximation corresponds to expanding exp(i(7oni) ~ 
1 + igont — .go„t^/2 + o{g^^). This results in Eq.(7) being 
rewritten as 



(exp(it ^ gonNn))T = 

n>0 

exp I I ^ gonNnt - t^/rl j , 

\ n>0 / 



(9) 



where the notation 



Td 



n>Q 



'N„) 



(10) 



for the dephasing rate is introduced. 

Eqs. (9) and (10) indicate that a collective mode of 
the confined atomic condensate should exhibit a dephas- 
ing induced by thermal fluctuations of the normal com- 
ponent. Below we suggest a model of dephasing of a col- 
lective mode of the confined condensate in the Thomas 
Fermi limit, and calculate the value of Td for the isotropic 
trap. 



III. ADIABATIC EFFECTIVE ACTION FOR THE 
LOW ENERGY COLLECTIVE MODES 

In order to obtain 1/rd, one should find the coeffi- 
cients gon and perform the summation in Eq.(lO). An 
explicit expression for gon can be derived from a many 
body Hamiltonian H taken in the standard form 






2m ^^2 

1=1.2,3 



2 2 

mcj, r , 



(11) 



where the Bose operators 4*, ^^ obey the usual Bose com- 
mutation rule; fj, is the chemical potential; tOi denote 
three frequencies characterizing the trapping potential; 
and the interaction constant Uo — 47r?i a/m is expressed 
in terms of the scattering length a and the atomic mass 
771. As usual, in the presence of the condensate, one uses 
the conventional representation 



-qj = $^ + *', 



fl2) 



where $c is a classical condensate wave function giving 
the condensate density 77c — l^cP and normalized to 
the number Nc of atoms in the condensate; $' stands 
for the excitation part. For the latter, the Bogolubov 
representation should be employed as 



^' = J2iUnan + V:ai), 



(13) 



where a„ destroys a quasiparticle on the level having the 
energy En, and (C/„, Vn) is an eigenvector of the linearized 
Bogolubov equations 



EnUn ^ H[Un + Uo^lVn, 
-E„Vn=H[Vn+Uo<^fUn 



(14) 



H[=Hi + 27io|*c 



/i. 



Expressed in terms of the quasiparticle operators ama\n 
the Hamiltonian (11) acquires the form 



H — 2_^ Ena\an + Hint , 



where the interaction part Hint contains the terms 

Hint = ^{gmnkamanak +H.C.) + 

mnk 



^ gmnkia'Lo'iakai, 

mnkl 



(15) 



(16) 



which could be identified with the terms ^ N^Nn in 
Eq.(5). Specifically, performing calculations in the first 
and in the second orders of the perturbation theory with 
respect to ,ff,-„(, one finds the coefhcients gkn in Eq.(5) 
as 



— An -^0\^\n |2r '^{En-Ena) 

gkn ^gknkn i ^ / ^ Igmnk] 



\En-EmY-El 



(17) 



En + En 



Ek 



where the coefficients gknkn and gmni can be expressed 
explicitly in terms of the Bogolubov amplitudes in 
Eq.(13) by means of substituting Eq.(12) in Eq.(ll) and 
selecting the terms (16). This yields 



gknkn = ^Jdr[{\Uk\' + \VkmUn\' + |KP) 

+ {U*kVkV:Un + c.c)] 



(18) 



and 



ffmnfe = ^0 / dr^,[Uk{U*^Un + V^^n + U*„yn) + 
Vk{U-:^Un + V;nVn+V::,Un)] 



(19) 



for real $c- 

We note that in the case of a continuum or quasi- 
continuum spectrum |18|| , the sum in Eq.(17) gives the 
main contribution to the imaginary part corresponding 
to the LD in the lowest order of the perturbation theory 
]l3|- p5| , p8[ . In the isotropic trap, this imaginary contri- 
bution is not significant ||l^ , which implies that the LD is 



suppressed. Therefore, the real value oi gkn (17) could be 
employed in Eq.(lO) for calculating l/r^. Unfortunately, 
such an approach leads to an expression l/r^ which for- 
mally diverges on low energies, despite the natural ex- 
pectation that the main contribution to the dephasing 
is produced by high energy excitations. Accordingly, the 
rate (10) acquires the incorrect T-dependence l/rd ~ VT 
for T > fj,. In fact, this divergence can be eliminated by 
a proper renormalization of the vertex in the low energy 
region, where an adequate description relies on the hydro- 
dynamical approach [ pl| . In order to solve this problem 
for the low energy collective modes, we employ a simple 
scaling procedure which yields a description of the effect 
of thermal dephasing in closed form at finite T in the 
limit of large N. 

First, we note that the dephasing effect is an adiabatic 
process when the high energy component follows the evo- 
lution of the low energy collective mode without dissipat- 
ing energy. As was discussed in Rcf. [g2|,|3[, the low en- 
ergy confined condensate modes at T = can be viewed 
as time dependent scaling r, -^ ri/bi of the coordinates r^ 
by some time dependent scaling variables bi = bi (t) . Fur- 
thermore, it has been shown [E2| that, if one ignores the 
kinetic energy term, such a scaling approach is exact for 
any given initial state of the many body wave function. 
As will be seen below, this implies that no dephasing of 
the low energy scaling modes should occur in such an 
approximation. The dephasing is induced by the kinetic 
energy terms, which are, however, small in the Thomas 
Fermi limit ||22| . Such a smallness implies that one can 
still consider the scaling variables bi as proper collective 
degrees of freedom, whose dynamics is modified in the 
presence of the kinetic terms. Below we will derive an 
adiabatic classical action for the 6-variables. In order to 
do this, we treat the ^-operator as a classical field by 



means of considering the 



it 



operators in Eq. 



as 



c-numbers. Then we employ the scaling ansatz 



*(n,t) 



V6IM3 



*(^,m), ^=fEr^^ (20) 



oscillate in time. The diagonal terms a*a„, which do 
not oscillate, should be retained and identified with the 
population factors. Let us denote the action obtained as 
a result of such a procedure as Sb = (S) . Then, we find 



where the notations 



R, - ( /"dr^Vf*), P, = ( /"dr^V,**V,*), 



(23) 



G^ (y / dr******) 



are introduced. Taking into account that these quantities 
do not depend on bi, one can vary Sb with respect to bi 
and obtain the classical equations of motion 



b, + Ujfbi 



G 



1 



P^]_ 

mRi 5^616263 mRi bf 



= 0, i = 1,2,3. 



(24) 



We note that these equations reproduce correctly the 
T ~ low energy spectrum obtained in Ref. [Q for the 
isotropic trap characterized by the relation uji = i02 = 
Lu^ = LOho- Indeed, in the case of zero temperature one 
should take the means (23) over the condensate state by 
setting a„ = aj^ = in Eq.(13), and, accordingly, taking 
^ = $c in Eq.(23). Then, the virial theorem yields 



miofR^"^ 



p{c) _ ^(,) ^ 



(25) 



with the superscript (c) indicating that the means in 
Eq.(23) are taken over the condensate state. Employ- 
ing this relation in Eqs.(24), one finds 



^il-^) 



rT^-d'Wi = o, (26) 

hbib2b:i ' ' W 



where the notation 



where t'{t) is some function of time determined in terms 
of the variables bi and their time derivatives bi. In what 
follows we assume that the scale invariant shape of ^ 
is given by Eqs.(12)-(14) obtained for bi — 1. In this 
manner we eliminate the non-adiabatic processes induced 
due to bi 7^ 0. Consequently, one can derive an effective 
classical action Sb = S[bi, bi] for the variables bi by means 
of performing the scaling transformation (20) in the full 
classical action 



S= dt{ 



dY-{^*-^ -h.c.)-H}. 



(21) 



Note that a substitution of the form (12), (13) is to be 
done in Eq.(ll) and, then, in Eq.(21). Then, the off- 
diagonal products a*^am and a„am {m ^ n) should be 
eliminated in the adiabatic approximation because these 






p. 



(c) 



(c) 



mufR) 



t(c) 



(27) 



is introduced. We note that Q determines the ratio of 
the averages of the kinetic energy to the harmonic poten- 
tial energy both taken for the i-direction. Accordingly, 
for the case of the isotropic trap considered in Ref. |2J] 

one can find Q"' = Ekin/Eho, where Ekin, Eho stand 
for the total kinetic and the total harmonic energies, re- 
spectively. Then, from the linearized version of Eqs.(26) 
we reproduce the frequencies of the quadrupolar mode 
luq — \/2uJho{'^+Ekin/ EhoY^"^ and of the breathing mode 
WA/ = w,io(5 - Ekin/EhoY/'^ derived in Ref. Q. 

Note also that for T 7^ 0, Eqs.(24) coincide with those 
derived in Refs. p3,E3[ in the Thomas-Fermi limit. This 
can be seen by means of setting the kinetic energy terms 



Pi = in Eqs.(24) and performing the scaling transfor- 
mation 



KijUi 



(28) 



with some k^ chosen in such a way as to make the so- 
lution 6i = 1 to be the equilibrium one. Furthermore, 
it can be seen that in the isotropic 2D case, when only 
two scaling variables 5i, 62 should be taken into account, 
the dependence on Pi can be eliminated by the scaling 
transformation (28), so that the frequency of the breath- 
ing mode does not depend on Pi. This implies that no 
dephasing of the 2D breathing mode should occur in ac- 
cordance with the result of exact calculations of Ref. [g5| . 
In order to simplify the following analysis of the T ^ 
case, let us consider a breathing mode in the 3D isotropic 
trap. Thus, we set b = bi = b2 = b^ and b = nb in 
the transformation (20) as well as in the action Si, (21). 
Then, varying SSb/Sb — 0, we obtain the nonlinear equa- 
tion describing the low energy adiabatic dynamics of the 
breathing mode in the presence of the normal component: 



&4 64 



0, C 



p 



mRujf^^K 



4. (29) 



where P = Pi + P2 + P3, R = Ri + R2 + R3 and k obeys 
the relation 



2 



p 



3G 



rriK'^R mn^R 



0. 



(30) 



An explicit form of ^ can be found in the limit P —t 0, 
which corresponds to neglecting the kinetic energy of the 
system if compared with the interaction energy. Setting 
P = in Eq.(30), one finds k and then, resorting back 
to Eq.(29), obtains 



(31) 



W^f (771^)1/5(304/5 



In order to express P, R, G in terms of the products 
a^ttn which, as was mentioned above, should be then 
identified with the population factors iV„ = a*^an of 
the quasiparticles in the second quantized picture, one 
should use the representations (12), (13), (23). How- 
ever, in the Thomas Fermi limit valid for large numbers 
Nc of atoms in the condensate, one may neglect cor- 
rections to P, G and Nc due to the excitations at not 
very large temperatures. Indeed, in the condensate state 
R r^ G ^ Ncri, where Vc stands for the condensate ra- 



.1/5 



26 



Therefore, 

.3/5 



R r^ G 



m 



7/5 



The 



dius Tc ^ Nc 

kinetic term P ~ Nc/r1 ^ Nc''' . High energy excitations 
produce changes 5R, SG, SNc, SP of P, G, Nc, P, re- 
spectively. Therefore, their relative contributions SR/ R, 
5G/G, SNc/Nc and SP/P to the value of ^ are very dif- 
ferent. Specifically, the first two terms contain an addi- 
tional smallness ~ Nc if compared with the last one. 
Therefore, in calculating ^ in Eq.(31), one should take 
into account only the contribution due to (5P, and take 



the values R, G determined for the condensate state by 
Eq.(23) for some mean value of Nc. 

The term SNc/Nc, arising due to the conservation of 
the total number of particles, can be neglected as well. 
Indeed, if the normal component fluctuates by some 
number 6N' , the conservation of the total number im- 
plies that 6Nc = —5N'. The kinetic term fluctuates as 
SP ^ 6N'. Therefore, the ratio SNc/Nc over SP/P con- 

— 2/5 

tains a smallness r^ Nc << 1. 

Finally, employing representations (12), (13) in (23) 



and taking the means so that {a] 



NnSmn, One finds 



C = ^^'^+E„5nA^n, 



m^RcUjl^ 



dr(C/;AC/„ + y„*AC/„), (32) 



Rc = fdrr^l^c 



where the virial relation (25) has been utilized, and ^'^'^^ 
denotes the contribution due to the kinetic energy of the 
condensate. In what follows we will neglect this term 
which does not depend on Nn- 

The solution b of Eq.(29) should be averaged over ^ 
represented by Eq.(32). Such an averaging can be un- 
derstood as a thermal ensemble averaging over possible 
Fock states of the quasiparticles. This interpretation 
closely resembles the case of destructive measurements 
pO| , when the initial conditions determined by Nn for 
each newly created condensate can vary from one con- 
densate to another. In the case of non-destructive mea- 
surements, the averaging of the solution b should be per- 
formed over a single many body state, which is a mixed 
state rather than a pure Fock state with respect to Nn- 
In accordance with the ergodic hypothesis such an av- 
eraging should give the same result as that obtained by 
means of thermal averaging as long as the number of 
quasiparticles is sufficiently large. In what follows we 
will not distinguish these two cases, and will employ the 
thermal averaging (...)t with respect to Nn- This averag- 
ing can be performed explicitly for the linearized solution 
of Eq.(29). Specifically, representing b = 1 + r; in Eq.(29) 
and keeping the terms linear in 77, one finds 



(^(0>T=%(e'"^°('-«''''*)T + c.c. 



, ,2/2 



(33) 



^u 



V^UJho 



2V5 



(e). 



where rjQ = const accounts for the initial condition 5(0) = 
1 -f 770 1 and the dephasing rate of the breathing mode is 



:^=^^oj^J29'nUl+Nnl 



TdM 



Nn^ 



40 



(34) 



exp(|-)-r 



with the coefficients 5„ given by Eq.(32). Performing 
similar calculations for the quadrupolar mode, we find 
the relation I/tv^q = \/5/2/TdM for the dephasing rate 
of the quadrupolar mode. Taking into account Eqs.(32), 
(34), we obtain an explicit expression for the rate I/t^a/ 
(34) in the WKB approximation |2^ (see Appendix A) 



as 



1 

TdM 



TmD{T/2ii), 



where the coefficient T m is 

2 



35 



M 



Tho 



a 



Tho 



-^ho, Tho 



mujfio 



(35) 



(36) 



and the universal dimensionless function -D(/3) is defined 
in Appendix A (see Fig.l), with the parameters Tc, /i 
given explicitly in (A2). 

In the limits /3 >> 1 (T >> 2^) and /3 << 1 
(T << 2/i) the function D{(3) can be found explicitly 
(see Eqs. (A28) and (A32), respectively). The current 
experimental situation is closer to the first case. It is con- 
venient to express T in units of the transition tempera- 
ture Tc of the Bose-Einstein condensation in the isotropic 
oscillator trap 



Tr — huj 



ho 



f ^ 



VC(3) 



1/3 



(37) 



where C(3) ~ 1.202; N is the total number of the trapped 
atoms (for T not very close to Tc we set Nc ~ N) . Then, 
we find 



1 

TdM 



35x/0.3 



15V5(C(3))5/6 



5/2 



a 



2/5 



N 



-17/30 



^ho 



(38) 



from Eqs.(35)-(37), (A28), (A29). Choosing the values 
T/Tc = 0.9, N = 2-103-104 andw/io = 27r-200s-\ rho = 
lO^^cm typical for the experiment [Q, we obtain the 
rate 1/TdM ~ 40s~^ — 20s~-'^. We note that these values of 
^/TdM are close to the damping rate observed in Ref. Q . 
However, for the chosen parameters, /3 = T/2fi « 1.3 
which is far from the requirement /? >> 1 insuring the 
validity of Eq.(38). Nevertheless, the above estimates re- 
main valid. Indeed, evaluating the complete expressions 
(A18), (A19) numerically (see Fig.l) changes these esti- 
mates by only about 20%. Specifically, the rate becomes 
« 50s~^ — 25s^^. For the lowest temperature achieved in 
the experiment |10| T « 0.4Tc, Eq.(38) becomes invalid 
because this temperature corresponds to /3 « 0.6 < 1. 



Accordingly, a numerical evaluation of D{(3) by means of 
Eqs. (A18), (A19) and, then, a substitution of the result 
into Eq.(35) yields the rate I/t^a/ = 8s~^ — 4s^^ which 
is also in the range obtained by the JILA group [[lO|. 




FIG. 1. The universal function D[j3): the solid line is the 
result of numerical calculation (A18),(A19); the dashed line 
is the approximation (A28). 

We emphasize that in the anisotropic trap employed by 
the JILA group the damping is most likely to be caused 
by the LD |12-nq,nq], and not by the mechanism dis- 
cussed above. The discussed mechanism in its pure form 
can be realised in the isotropic trap only. Therefore, a 
correspondence between the rate calculated above for the 
isotropic trap and that measured by the JILA group |l^] 
for the anisotropic trap indicates that, while decreasing 
a degree of the trap anisotropy, the damping rate should 
practically stay unchanged despite the fact that the na- 
ture of the damping changes. 

In the case T << 2fj, one could use an explicit form 
(A32) for D{P) and, correspondingly, find an explicit T- 
dependence of the rate (35). However, in this case the 
rate becomes so small, that the mechanism of the quan- 
tum self-dephasing MM comes into play. 

It is interesting to investigate the dependence of the 
dephasing rate on the amplitude of the oscillations. It is 
worth noting that in the case T = such a dependence 
is very pronounced MM- In contrast, as will be seen 
below, the amplitude dependence at finite T is weak. In- 
deed, this dependence is due to the nonlinearity of the 
term ~ ^ in Eq.(29). In the lowest order with respect to 
the initial value r/o in Eq.(33), this dependence can be ob- 
tained by expanding Eq.(29) up to the terms r^ rf and r/^ 
and finding the correction to the frequency of the lowest 
harmonic in the order ~ I^P- Performing straightfor- 
ward calulations (see Appendix B) and then averaging 
over the ensemble, we obtain 



TdM 



TdM{Al) 



llA.l' 



(39) 



the ratio of the rate l/rdM(^i) determined in the second 
order with respect to the amplitude Ai — 2r]o of the 
collective mode to the rate I/tclm in the 0-th order given 
by Eqs.(34), (35). As one can see, the rate demonstrates 
a slow decrease as a function of the amplitude Ai « 1. 
Here we have shown that the collective excitations 
of the confined Bose-Einstein condensate should demon- 
strate a dephasing caused by thermal fluctuations of the 
normal component. In the following section we will dis- 
cuss how this dephasing effect can be distinguished from 
irreversible dissipation experimentally. 



IV. THE ECHO EFFECT IN A CONFINED 
BOSE-EINSTEIN CONDENSATE 

The reversible nature of the damping can be tested in 
an echo experiment similar to the spin echo, the pho- 
ton echo and the phonon echo effects (see in Ref. [pO|). 
The nature of this effect can be briefly outlined as fol- 
lows poi . A short external pulse imposed on the system 
at the time i = excites a collective mode. The collec- 
tive mode amplitude decays due to dephasing as well as 
due to irreversible dissipation. Both processes are char- 
acterized by their typical rates IJTd and 7, respectively. 
The second pulse imposed at the time t — t partly re- 
verses in time the evolution of the system initiated by 
the first pulse. This implies a partial revival of the de- 
phased amplitude at the time t « 2r. We note that the 
occurrence of the echo is a general property of the system 
where irreversible damping is weaker than the dephasing. 
Thus, a necessary condition for observing a distinct echo 
is Td < 1/7, Td < r < 1/7 . 

Specific features of the echo depend on the details of 
the system. The time profiles of the responses, as Eq.(33) 
indicates, should be gaussian in the case of the thermal 
dephasing discussed above. In the case of the LD these 
responses should be characterized by exponential relax- 
ation. Presently available experimental data |jl3,pi do 
not allow the distinguishing of the gaussian type damp- 
ing from the exponential one [GSJ. In the next paper we 
will analyze the echo in the anisotropic confined conden- 
sate, where the main cause of the damping is the LD. 
Below we will study the situation in the isotropic trap, 
where the dephasing is caused by the thermal mechanism 
described in Sec. III. 

A relevant description for the case under considera- 
tion relies on Eq.(29) modified to incorporate the exter- 
nal drive as well as some possible irreversible dissipation. 
As was discussed in Ref. Q, the external drive (5u;^(i), 
which changes the curvature of the trapping potential, 
should be included in the linear part of the equation for 
the scaling variable b. Accordingly, Eq.(29) is rewritten 
as 



b+[ul + Su;Ht)]b- 



64 



27& + iuj_ 



1 



ho ^4 



= 0. (40) 



For £_= 7 — 0, one obtains the equation derived in Refs. 
p2| , P3| for the case T = 0. The term '--^ 7 describes 
the irreversible dissipation at T 7^ 0. The term ^ ^ 
introduced already in Eq.(29), with ^ given by Eq.(32), 
accounts for the dephasing effect discussed above. 

The time dependent part 5uP'{t) of the frequency 
should be driven so as to be in resonance with the col- 
lective mode, that is in the form 

6u\t) - -/(i) exp(ic^oO - r (t)exp(-iwoi), (41) 

where wq = V^^ho and f{t) stands for the complex am- 
plitude of the external drive. In order to avoid exciting 
other modes of the system, this amplitude should be con- 
sidered as a slow envelope of the resonant drive with a 



typical time t/ >> 



The echo, then, can be pro- 



duced by making J{t) reach a maximum at i = and 
then become zero untill the time t = t, when f{t) peaks 
again. For sake of simplicity, we will ignore other modes 
and will analyze the simplest situation when the external 
drive produces two (5-pulses 



Suj\t) = -fi5it)-f2S{t-T) 



(42) 



at i = and t — t having amplitudes /i, /2, respectively. 
For the case of small amplitudes /i, /2 of the drive, 
one should look for an evolution of the small perturba- 
tion around the equlibrium value 6=1. We note that, 
in contrast with the conventional situation J20| , the echo 
response in our model does not require non-linearity of 
the dynamical equation. This is due to the fact that 
the external drive plays a two-fold role. Specifically, 
on one hand, it gives rise to an effective external force 
—5uj^{t) and, on the other hand, it excites the system 
parametrically. Indeed, linearizing Eq.(40) by the sub- 
stitute 6 = 1 + ?7, with 77 << 1, one obtains 

V + kL(5 -0 + 3ojHt)]7^ + 277? = -Suj^t), (43) 

where the higher order terms of 77 are neglected. 

We assume that initially ai t = —00 the mode was 
not excited (77(— 00) = 77(— 00) = 0). Then taking into 
account Eq.(42), one finds form Eq.(43) 



77(0)^0, 77(0) = /i 



(44) 



after the first pulse. The second pulse at t = t results in 
a jump of 77 so that 



77(r + £)=77(T-£) + /2(l+77(T)), 

■nir + e) =7/(T-e) = 77(t), 



(45) 



where s ~> +0. 

We are looking for a solution at f > t. It has the form 



77(t) = ^e(*'3"'^)(*^^' + yl*e(~''^^'^'(*~'^) 

l^ ' (46) 

lO'^' 



LOoi5'0'^'^V5LOo{l-h, 



where we have taken into account that 7 << uJho and 
^ << ujfio- An explicit expression for the coefficient A 



can be obtained if one employs the conditions (44), (45). 
Finally, we find the solution (46) for i > r expressed as 



v{t)^Ml + 43)^^"'-'^' 






Veit) 



C.C. 



where 



.X _ /2/1 »Q(t-2r)-7t 

represents the echo occurring at the time moment t = 
Indeed, after thermal averaging over Nn, one finds 



(47) 



(48) 



2r. 



{v{t)h 



A 

UJq' 



[sin(a)t) 



A 

2uj< 



cos(a;t)]e 



-tt-tVr^, 



/2/: 



+ ^ sin(w(i - T))e-^(*-")-(*-")'/"'A^ (49) 

UJQ 



10^1 



^ cos(w(t ~ 2T))e-^*-(*-2-)'/^" 



The echo amplitude can be defined as Ae — {ri{2T))T in 
the case when the last term in Eq.(49) dominates the 
sum. In this case we find 



. = J2A ^.2,. 



lo^L 



(50) 



In deriving Eqs.(49), (50) in the limit under considera- 
tion, we have made the replacement Q — \/^uJho every- 
where in Eqs.(47), (48) except in the exponents, where Q 
linearized in S, and given by Eq.(46) has been employed. 
Thus, we obtained the echo effect in the linear approxi- 
mation. 

We have also analyzed the non-linear echo problem for 
Eq. (40) numerically. This equation was solved for a given 
value of ^, and then the final solution was averaged over 
the values of f distributed in accordance with the gaus- 
sian G{£) — ^/•v/7rexp(— <^^/0^), where 9 determines the 
effective width of the distribution in such a way that the 
averaging of the linearized solution reproduces the result 
(33), (34). Specifically, we set 9 = V^/uihoTdM- The 
results of the calculations are shown in Fig. 2. In the case 
a) the amplitude of the second pulse is too small to make 
the echo observable. In the case b) the second ampli- 
tude is 5 times stronger so that the echo is distinct. In 
the case c), while the second pulse amplitude /2 is the 
same as in the cases a) and b) , the amplitude of the first 
pulse /i is two times larger than that in the cases a),b) 
(note the different scale of the vertical axis in the case 
c)). As one can see, the echo in this case merges with the 
tail of the second pulse, which creates an impression that 
the decay time of the second pulse increases by several 
times. In order to produce the echo in the case of the 
large amplitudes /i, /2, the time separation r between 
the pulses should be increased. However, in this case the 



irreversible dissipation may strongly suppress the echo in 
accordance with Eq.(50). 




FIG. 2. The echo effect in a confined Bose-Einstein con- 
densate: the solution for the breathing mode h{t) given by 
Eqs.(40), (42) ( 7 = 0.01aj,,o, 6^ = 0.02, r = 80a;^„'). 
The cases a), b) and c) are different in amplitudes /i, /2: 
a) /i = O.^ujho, /2 = 0.1c<;;io; b) /i = /2 = O.^cJho; c) 

/l = luJho, /2 = 0.5iJho- 

The echo effect analyzed above is a classical mechan- 
ics effect. Below certain temperature Tq, the rate of 
the quantum dephasing |g,|[ should become faster than 
the damping induced by the normal component. Accord- 
ingly, the classical treatment employed above becomes 
no longer valid. The problem should be reformulated in 
terms of the quantum dynamics of the variable 77 in a 
sense of the approach |^, with the external drive (42) 
taken into account. It can be shown that the echo still 
exists at t = 2r. Therefore, the spontaneous quantum re- 
vival, determined by the interaction constant and thereby 
occuring at very long times H , can be induced to occur 
at much shorter times comparable with the time of the 
quantum collapse H . This problem will be considered in 
a separate publication. 



V. DISCUSSION 

We have suggested a mechanism for the apparent 
damping of a Bose-Einstein condensate confined in the 
isotropic oscillator trapping potential. This damping is 
a reversible dephasing of the collective modes caused by 
thermal fluctuations of the population factors of the nor- 
mal component. The calculation of the dephasing rate 
gives a value which is comparable with the experimen- 
tally observed rate of the damping of the low energy col- 
lective modes in the atomic traps. 

This mechanism of dephasing relies on the ensemble 
averaging of the collective mode over the initial popula- 
tion of the normal component. Thus, an assumption is 



made that for any given initial distribution of the popula- 
tion factors of the "hot" quasiparticles, this distribution 
does not relax to equilibrium during the time of the de- 
phasing Td- Accordingly, processes of relaxation due to 
the LD or collisions may suppress the discussed mecha- 
nism, if their relaxation times are comparable with r^. As 
long as the coUisional damping introducing irreversibil- 
ity is unlikely to be relevant for such small temperatures 
and densities, the LD is the only competing mechanism. 
However, the LD is expected to be significant for sub- 
stantially anisotropic traps only [Q. Therefore, in the 
traps characterized by small anisotropy our mechanism 
should dominate. 

Both mechanisms of damping - Landau damping and 
that considered above - are reversible in nature, and 
therefore the evolution of the system can be partly re- 
versed in time. We suggest testing this in the atomic 
traps by employing the echo effect. As our analytical 
and numerical calculations for the breathing mode in the 
isotropic trap indicate, the echo amplitude as well as its 
position depend on the parameters of the external drive 
which can be varied over a wide range. 
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APPENDIX A: WKB CALCULATION OF THE 
DEPHASING RATE 

The WKB calculation of the dephasing rate presented 
here is essentially based on the results of Rcf . |2^ . Em- 
ploying Eqs.(14) as well as the normalization condition 
! dv{\Un? - IKP) = 1 in Eq.(32), one finds 



mRcUjl^ 



{En- 



(Al) 



{KVnU;, + c.c.)]}, 

where the notation K = Uo$^ is introduced, and for 
the condensate wave function $c ~ v^ we employ the 
Thomas Fermi solution lEQl 



nc='^^{rl-r'^)Q{rc-r), 



1/5 



re = rh, 



V »■'- ) ' 



(A2) 



were Q{z) is the step function; Uo and Vho are defined in 
Eq.(ll) and Eq.(36), respectively. Accordingly, one finds 
the value of Re in Eq.(32) as 



Re 



SSa'^L ' 



(A3) 



States in the isotropic trap can be classified in terms 
of the angular momentum L, its z-component Lz and the 
radial quantum number n^. Thus, the index in (Al) as 
well as in the sum (34) should be understood as consist- 
ing of these three quantum numbers. This implies that 
the summation ^^^ ... in Eq.(34) runs over three quan- 
tum numbers n = {nj.,L,Lz). Because of the spherical 
symmetry, the summation over Lz can be performed triv- 



ially, which gives ^^ 



En.EL(2i + l).... Aswillbe 



seen below, the large values Ur » I, L » 1 dominate 
this sum. Therefore, we replace the summation by the 
integration over n^, L 



drir 



La 



dL2L. 



(A4) 



where the upper limit Lq is to be determined, and we 
made the replacement 2L -t- 1 « 2L. It is convenient to 
change the variable n^ to E by employing the quantiza- 
tion condition IETI 



1 1 r^ 

nr + - = ^ / drpr, 

2 TTh ,L, 



Pr = ^2m{^/WT\K?~Ueftir)). 



(A5) 



where 



Ueff{r) = -7nujiy + 



n^{L + l/2f 
2mr^ 



2\K\-^i (A6) 



denotes the effective WKB potential ^j , and the turning 
points ri , r2 obey the equation pr = or 



VE'^ + K^-Ueff{r) =0. 
Then the integral (A4) acquires the form 



2 f°° f^" r^ dr 

4 / dE I dLL I —. 



"'0 



(A7) 



(A8) 



where Eq.(A5) has been employed, and Vr stands for the 
WKB radial velocity H]. 

Before we proceed, it is convenient to employ dimen- 
sionless variables of length, energy and angular momen- 
tum as 



re riLOho ri ri 



(A9) 



respectively. Note that in these units the condensate ra- 
dius Tc equals 1, and the chemical potential /i and the 
quantity K become 



respectively. Accordingly, Eq. ( A7) yields two sets of solu- 
tions for the dimensionless turning points a;i,2 = tx^iJTc 



Xi 



y/y^, X2 = v^> J < ^2e 



(All) 



where 



2/0 



j'^(l + Vl + 4e2 + 2J2) 
4e2 + 2 J2 ^ 



(A12) 



and 



6± = ^+^±V(^+>-^^ 



Xi = \/&-' 2^2 = 



2- 



(A13) 



/2e < J < e + i e > 



As has been discussed in Ref. [^, the solutions (All) 
and (A13) correspond to the case when the classically al- 
lowed region extends into the condensate and to the case 
when it is totally outside the condensate, respectively. 

The U, y-amplitudes inside the classically allowed re- 
gion are M 



U^^2(^iS^ + S^)^Y, 



2r\ 



3/2 



X^Vx 



L,L,, 



F = ^^(5+-5_)^y,,,. 



(A14) 



2r: 



3/2 



harmonic: the normalization constant is 



Here S± = a/-\/1 + (k/e)^ ± (fc/e); Yl^l^ is the spherical 

(A15) 



and th e dimensionles s radial velocity 
Vx = yjmr'l^/hujhorlvr is given by 



Vx ^ \ 2 



fc2 



(y/JTk^-k-^), 



(A16) 



Xl < X < 1, 



inside the condensate and by 



v/2e+l-a;2- j2/x2, 
1 < X < 0:2 



(A17) 



outside the condensate. In calculation of Co(e, J) and in 
what follows, we replace sin2 ip by 1/2 because the WKB 
phase (j) [pTl varies rapidly inside the classically allowed 



region. The integrals outside this region are exponen- 
tially small, and we neglect them. 

Substituting (A14) into (Al) and employing the units 
(A9) in Eq.(34), we find the expressions (35), (36) where 
the dimensionless function D{P) is defined as 



Dm 

with the notation 
Pie) 

1 



de- 



e/(3 



-,1/2 



rP(£) 



(e'//3 - 1) 



2 Jx^ Vx 



(A18) 



(A19) 



2 -'^ I r,7 ^ /1 I , ^ n 2 



(-x'-- + 2k)xl + - + -]} 



introduced, and k determined in Eq.(AlO). The value 
of the limit Jo(e) can be found from Eqs.(All), (A13). 
Specifically, for e < 1/2, only the case (All) can be real- 
ized. This implies that 



Me) 



v^e, €<-. 



For e > 1/2, Eq.(A13) yields 
Me)^e+-, 



e > 



(A20) 



(A21) 



Consequently, the integral (A19) can be expressed as 
p{e) = pi{e)+ p2{e), 

Pi{e) = j ''dJJ^^^^[Ini(e,J) + In2(l,6,J)]2, 

P2(e) = 

1 /•c+1/2 fi2( T\ 

e(._i) / dJJ^l^[In2(v/5r,.,J)]2, 

2 Jv2i 2 



(A22) 



where the notations 



Ini(e,J)=/ - 



In2(a,e, J) 



^^ dx eVe2 + fc2 



xi < 1; 



c fc + Ve2 + A:2 ' 

^+ dx , 1 1 2x 

— (e+ O ^ O^ )' " - 1' 
Ut 2 2 



(A23) 



are introduced; Co(e, J) is given by Eq.(Al5), and Vx 
is determined by Eqs.(Al6), (A17). Note that here we 
have employed the explicit expressions (All)- (A13) for 
the turning points. 

The value of the normalization constant Cq can be 
found explicitly [M. The integrals (A23) can also be 
calculated explicitly. We find 
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for e > 1/2, ^/2l < J <e + 1/2, and 
Cg (e, >/) Z' 2e arccos ai 



(A24) 



%/2i 



,/2 



arccos a2 



Ini(e,J) 



e arccos ai 1 



V2e2 + J2 4 
1 



V2e- J2^ 



2\/2 



In as, 



(A25) 



In2(l,e, J) = -(e+ -)arccosQ;2- 
for J < \/2e, where we have introduced the notations 



ai 



a2 



'2e2 + J2 _ £ + eVr+4e2 + 2j2 
2eVl + 4e2 + 2 J2 



i-6+J(6+i)2-j2 



\ 2J(e+i)2-j2 



(A26) 



a-i = -^[\/l + 2e - Vl + 4e2 + 2^2- 

v2 



1 + 2e + Vl + 4e2 + 2 J2] 



(l + 4e2 + 2j2)-i/4. 

We note that in the formal hmit (3 » 1, the func- 
tion D{j3) given by (A18) can be found exphcitly. In- 
deed, in this case the main contribution to (A19) comes 
from e >> 1. This imphes that only the term p2{f) in 
Eq.(A22) should be taken into account because it gives 
the highest power of e as /92(e) ^ e*. As simple analy- 
sis of (A22) shows, the term pi(e) ^ e^ . Thus, taking 
p(e) w p2(e) and combining Eqs. (A24), (A22), (A25) 
and (AI2), we find 



.(e) ^T^(.^ 



-fQ{e- 



16 



(A27) 



for e >> 1. Substituting this into (A18) and taking the 
limit (3 » 1, we obtain 



D{P) 



5/2 



(A28) 



which yields Eq.(38). This expression is shown in Fig.l 
by the dashed line. As one can see, in the range of (3 of 



the order of 1 the approximation (A28) underestimates 
the rate by approximately 20%. 

We note that actual values oi j3 — T/2/i are far from 
being (3 >> 1. Indeed, employing Eqs. (A2) and (37), 
we find 



''-i-(Tg^)""r"=(3)«- 



1/15 






(A29) 



which yields values /3 « 1 for the experiment [Q for 
T Ri Tc- Therefore, for these values the function D{P) 
should be found numerically (see the solid line in Fig.l). 
In the opposite limit /3 — + 0, which corresponds to 
T << 2fi or large N, the contribution to D{l3) due to 
P2(e) becomes exponentially small. Thus, the term ^ 
Pi{e) (A22) dominates in (A18). Taking into account 
that the effective values of e ~ /?, one may perform an 
expansion in terms of the small parameter e in Eqs. (A26) , 
(A25) and obtain 



p{e)^pi{e)=poie^^^, e ^ 0, 



where the notation 



Poi 



7r/2 



1 f"'- . [X 

—^ I ax sm X - 

V2J0 



[x -I- -j^ sin 2x(7 -1-11 cos2)] 



2M2 



(A30) 



(A31) 



X + h sin 2x 



has been introduced. A numerical evaluation of this in- 
tegral gives poi ~ 1.5. This yields for (AI8) 



Dip) « Dop'/\ 



Do 



poij^dx-. 



,7/2 



1/2 



(A32) 



4.4 



in the limit (3 << 1. 



APPENDIX B: CALCULATION OF THE 

AMPLITUDE DEPENDENCE OF THE 

DEPHASING RATE 

Expanding Eq.(29) up to the third order with respect 
to r/, one gets 



fj + LoliT] - aMff + PmV^ = 0, 
where the notations are 

< = (5 - ex, ttM = 10(1 - lex, 

(3m = 20(1 - IX. 



(Bl) 



(B2) 



The solution of (Bl) up to the second order with respect 
to 770 has a form 



V=^\Vo\' + {Voe'-' + c.c.y 
^(%^e»^- + c.c.). 



(B3) 



3uj 



M 
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where the effective frequency lo in the same order is 



W = LOM +w'|?7oP, 



ba\„ 3/3^/ 



'M 



3^Af 



2ujM 



(B4) 



Now employing (B2) in (B4) and performing the thermal 
averaging of (Bl) over ^ in the limit ^ << 1, we obtain 

{t^)j, = T^^e''^-')Tt~i'/rlM(M) + ^,,^^ (B5) 

where the constant as well as the second harmonic have 
been omitted; the dephasing rate l/TrfA/(^i) as a func- 
tion of the amplitude Ai = 2770 of the first harmonic is 
given by Eq.(39). 
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